van der Waals interaction in iron-chalcogenide superconductors 
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We demonstrate that the inclusion of van der Waals dispersive interaction sensibly improves 
the prediction of lattice constants by density functional theory in iron-chalcogenides (FeCh) super- 
conductor compounds, namely FeSe and FeTe. We show how generalized gradient approximation 
(GGA) for the exchange correlation potential overestimates the out-of-plane lattice constants in 
both compounds when compared with experiments. In addition, GGA predicts a too weak bonding 
between the neutral FeCh layers, with a sensible underestimation of the bulk modulus, van der 
Waals corrected simulations completely solve both problems, reconciling theoretical results with 
experiments. These findings must be considered when dealing with theoretical predictions in FeCh 
compounds. 

PACS numbers: 74.70.Xa, 71.15.Nc, 61. 50. Ah, 62.20.-x 
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I. INTRODUCTION 

The discovery of high-temperature superconductivity 
in iron pnictides (FePn)ii^ raised a strong interest in 
searching for new superconducting materials which con- 
tain Fe and share common structural features with FcPn 
superconductors. For example, superconductivity was 
found in the iron chalcogenide (FeCh) FeSe and its al- 
loys when Se is partially substituted by Te<^ The struc- 
ture of FeCh is characterized by stack of neutral layers 
tetrahedrally coordinating Fe ions with chalcogens (Ch), 
similar in structure with respect to FeAs layer in FePn's 
which, on the contrary, are negatively charged due to 
the presence of intercalates and oxides layerS)^ Despite 
lower critical temperature, the Fei+ySexTei-x family has 
a simpler crystal structure compared to those of FePn. 
For these reasons, the FeCh alloys can be considered as 
useful prototype systems to investigate the fundamental 
aspects of structural, electronic, magnetic and supercon- 
ducting properties in Fe-based superconductors. 

Moreover, FeCh alloys show exceptional physical prop- 
erties originating from the competing magnetic and su- 
perconducting orders. In particular, FeSe is supercon- 
ductor with a transition temperature of Tc ^ 8 K at 
ambient pressure;^ which grows up to 37 K when the 
pressure reaches 9 GPa;^ suggesting that the lattice 
plays a fundamental role in the superconducting tran- 
sition. On the contrary, the parent FeTe crystal is not 
supcrconductorfii^ However, it was shown that 90 nm 
thick films under tensile strain become superconductors 
with onset temperature at 13 K, confirming a sensible 
role of the lattice »^ 

First principles density functional theory (DFT) is con- 
sidered a fundamental theoretical and computational tool 
to investigate the structural and electronic properties of 
the normal state of Fe-based superconductors.— In- 
deed, it predicted the experimentally confirmed topology 
of the Fermi surfaces (FSs), the magnetic phases and the 
structural distortion observed at low temperature. 

In the early period of research, a strong interest was de- 
voted to understand the predicting power of DFT in both 



local density (LDA) and generalized gradient approxi- 
mations (GGA)fi^ showing the importance of magnetic 
correlation in the prediction of lattice constants. It was 
found that geometry optimizations performed consider- 
ing a static Fe magnetic moment (although much higher 
than measured one) generally mimic the magnetic fluctu- 
ations and correlations)^ thus predicting crystal lattice 
constants in acceptable agreement with experiments^^ 

However, a notable exception exists: ab-initio DFT in 
both LDA and GGA approximation fails in predicting 
the lattice constants of FeCh crystal structures, in par- 
ticular of the out-of-plane lattice constant Indeed, at 
present, most calculations have focused on the study of 
electronic properties considering the experimental mea- 
sured lattice parameters with the only optimized quan- 
tity being the Ch height (hch) with respect to the Fe- 
atoms plane4 '^°i^^ However, it must be emphasized that, 
due to the very simple layered structure, the equilibrium 
volume and hch, critically affect the electronic structure. 
Then, in order to predict the details of electronic prop- 
erties of FeSe, FeTe and their alloys from first principles 
in different physical conditions (alloys, pressure effect, 
surfaces etc.), it results fundamental to properly solve 
the theoretical problems affecting the out-of-plane inter- 
action between layers. 

Although widely recognized)^ the role of dispersive 
van der Waals (vdW) interactions between neutral FeSe 
and FeTe layers was not properly investigated. 

It is well known that a general drawback of all common 
exchange and correlation functionals is that they do not 
properly describe long-range electronic correlations, as 
the vdW interactionfii"— In fact, computational inves- 
tigations using DFT are not simply interpretable when 
studying systems in which vdW dispersion plays a crucial 
role due to non-local correlation effects. 

Recently, a large effort was devoted to take into ac- 
count vdW interactions and very interesting works were 
dedicated to point out the state-of-the-art in advanced 
materials.— For example, Hyldgaard, in an extensive 
review^ evaluated and established the limits and ranges 
of applicability of many different computational ap- 
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proachcs developed to account for vdW on a large family 
of materials ranging from insulators, semiconductors and 
metals. 

In the present work, we show how first principles DFT 
successfully describe crystal structure of FeSe and FeTe, 
with unprecedented agreement compared with experi- 
ments, once corrected to include the non-local vdW inter- 
action, improving the calculation of out-of-plane lattice 
constant, interlayer binding energy and bulk modulus. 
In addition, we show the effect of corrected lattice pa- 
rameters in the electronic properties (band structure and 
Fermi surface) in both FeSe and FeTe. 



II. COMPUTATIONAL DETAILS 

The calculations were performed using the Vienna Ab- 
initio Simulation Package (VASP)2ii^ within the gen- 
eralized gradient approximation (GGA)<^ The Perdew, 
Burke and Ernzerhof (PBE)24 functional was used to cal- 
culate the exchange-correlation potential. The GGA ap- 
proximation correctly predicts the ground state magnetic 
phase for FeSei^ and FeTci^ Winiarski et. al.^ showed 
that LDA estimates more precisely FeSe lattice constants 
than GGA, however the calculations were performed in 
the PbO-type tetragonal non-magnetic phase, neglecting 
fundamental structural distortions and magnetic effects. 
However, the relative success of the LDA in predicting 
the high temperature lattice parameters is merely due to 
an accidental cancellation of errors between the correla- 
tions and exchange energies <^ 

In this paper, we used projected augmented-wave 
(PAW) pseudopotentials^ for all the atomic species in- 
volved and in order to achieve a satisfactory degree of 
convergence the Sd^ states of Fe, 4s^ of Se 
and 5s^ states of Te were treated as valence electrons 
with an energy cutoff up to 550 eV. Integrations over the 
Brillouin Zone (BZ) was performed considering different 
uniform Monkhorst and Pack grids^^ depending on lat- 
tices: 13x13x9 and 14x7x9 for magnetic coUinear stripe 
(AFMl) FeSe {a\/2 x b\/2 x c unit cell), and magnetic 
bicollinear double stripe (AFM2) FeTe (a x 26 x c crystal 
unit cell), respectively. For the tetragonal (a x a x c) 
paramagnetic (PM) phase, which contains two Fe and 
two Ch atoms, we used 20x20x15 and 15x15x9 k-grid 
for FeSe and FeTe, respectively. 

The vdW interaction is considered using the DFT-D2 
Grimme's semi-empirical force-field correctio n^"'^^ and 
the so-called vdW-optB86b functional as implemented 
in the VASP codci^ The two functionals were chosen 
for their simplicity (DFT-D2) and high accuracy (vdW- 
optB86b). 

Due to the high accuracy required in the calculations, 
we previously checked the pseudopotential quality with 
all-electron full potential linear augmented plane-wave 
method in the FLAIR implementation<^i2^ In Fig. [1] we 
show the c-axis relaxation, fixing the in-plane lattice pa- 
rameter to experimental values (see below and Tab. |T]for 



all details and references) on AFMl FeSe for VASP and 
FLAIR simulations. The pseudo-potential energy curve 
nicely agrees with the all-electron one in a wide range of 
c lattice constant, giving equilibrium c of 6.30 and 6.25, 
respectively. This is a fundamental consistence check 
due to the already discussed issues related to the com- 
parison between all-electron within the well converged 
pseudopotentialsji^ 
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FIG. 1: The pseudopotential (open circles) and all-electron 
(full triangles) energy curves calculated for FeSe AFMl phase 
as a function of c lattice constant. The in- plane a and b lattice 
constants are kept fixed at the experimental values.— The 
solid black line is a guide for the eyes. 



III. RESULTS AND DISCUSSION 

Having tested the accuracy of pseudopotentials, we cal- 
culated the energy-volume phase diagrams for FeSe and 
FeTe in both PM and AFM configurations with and with- 
out the vdW correction. The results (with available ex- 
perimental measurements) are presented in Tab. |T] and 
Fig. [2] and [31 

Since the PM calculations neglect magnetic interac- 
tions, fundamental to reproduce experiments, we will 
mainly focus our attention on AFM phases. 

In both FeSe and FeTe the GGA curves show a very 
weak interaction between the layers, predicting a too 
large lattice c constant when compared with experiments. 
As we can see, for FeSe (FeTe) the GGA gives a ~15% 
(10%) deviation from experiments for the out-of-plane c 
lattice constant, while the in-plane a and h parameters 
are in good agreement with a deviation lower than 1% 
(4%). 

On the contrary, the interlayer bonding energy, cor- 
rected with the semi-empirical DFT-D2 vdW dispersion 
potential, shifts the minimum of the total energy indicat- 
ing an increased interaction between FeCh layers, shrink- 
ing all the lattice parameters with respect to GGA. In 
FeSe and FeTe a and h remain consistent with experimen- 
tal values in a range of ^0%-5%, depending on material. 
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TABLE I: Structural parameters and magnetic moments for 
FeSe and FeTe. Available experiments are also reported. 



EXP 


T (K) 


a (A) 


b (A) 


c(A) 


hch (A) 


a (deg) 


FeSe2£ 


7 


3.7646 


3.7540 


5.479 


1.4621 


90.00 


FeTc^^ 


2 


3.8312 


3.7830 


6.264 


1.7540 


89.17 


GGA 




a (A) 


b (A) 


c (A) 


hch (A) 




FeSe 


PM 


3.68 


3.68 


6.26 


1.39 


90.00 


FeTe 


3.81 


3.81 


6.52 


1.59 


90.00 


FeSe 


AFM 


3.75 


3.71 


6.32 


1.45 


90.00 


FeTe 


3.87 


3.63 


6.90 


1.78 


86.60 


DFT-D2 




a (A) 


fe(A) 


c(A) 


hch (A) 


a (deg) 


FeSe 


PM 


3.64 


3.64 


5.42 


1.40 


90.00 


FeTe 


3.77 


3.77 


6.03 


1.59 


90.00 


FeSe 


AFM 


3.67 


3.61 


5.53 


1.46 


90.00 


FeTe 


3.81 


3.61 


6.42 


1.77 


88.53 



while along the c-axis, there is a sensible improving: the 
theoretical value is corrected within 1-3% with respect 
to experiments. Moreover, it is very interesting to note 
that, for FeTe AFM phase, the DFT-D2 improve also the 
monoclinic a angle. 

We have fitted the curves with a Birch-Murnaghan 
state equation,— and compare the equilibrium volumes 
and the bulk moduli (shown in Tab. with available 
experiments. 
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FIG. 2: FeSe energy-volume curves for GGA (open circles) 
and vdW DFT-D2 (full squares) in the PM (upper panel) and 
AFMl (lower panel) phases. Volume and energies refer to the 
conventional unit cell. Solid lines are the fitted Murnaghan 



TABLE IL Equilibrium volumes and bulk moduli of the con- 
ventional cell calculated for paramagnetic (PM) and antifer- 
romagnetic (AFM: AFMl and AFM2 for FeSe and FeTe, re- 
spectively) phases with and without vdW correction. 



EXP 


T (K) 


v., ik') 


Bo (GPa) 


FeSe^a 


50 


77.56 


33 


FeTe^ 


300 


91.98 


36 


GGA 




(A^) 


Bo (GPa) 


FeSe 


PM 


84.95 


5.28 


FeTe 


94.73 


9.81 


FeSe 
FeTe 


AFM 


87.82 
96.65 


3.37 
9.71 


DFT-D2 




v., (A^) 


Bo (GPa) 


FeSe 


PM 


71.73 


37.57 


FeTe 


85.53 


36.46 


FeSe 
FeTe 


AFM 


73.18 
88.08 


34.67 
38.99 



The most evident result is the striking disagreement 
between GGA and experimental bulk moduli for both 
FeSe and FeTe. In fact, as evident from Fig. [5] and [31 the 
out-of-planc interaction is too weak resulting in a very 



low Bq. This behaviour is dependent on the magnetic 
phase considered. Interestingly, we observe that the vdW 
correction completely changes the physics and chemistry 
of the out-of-plane interaction thus resulting in a much 
better agreement with experiments. We note that, even 
in this case, both PM and AFM phases are corrected in 
the same way. The overall satisfactory agreement indi- 
cates that the vdW interlayer interaction is fundamen- 
tal to correctly reproduce the properties of FcCh com- 
pounds. To the best of our knowledge, the bulk moduli 
of both FcSc and FeTe have never been predicted so far. 

As evident from our results, and as already well dis- 
cussed in review articles (see, for example, Ref. [T7llT8ii2g ) . 
the main effect of the vdW interaction is the c lattice 
constant reduction and the consequent increase of Bq, an 
effect common to other layered materials. 

In recent reviews 3i it was shown as different ap- 
proaches to include the vdW interaction can lead to dif- 
ferent calculated lattice constants. The DFT-D2 is a 
semi-empirical method, very efficient, but relies on the 
optimization of four semi-empirical parameters,— previ- 
ously fitted on different classes of materials.— In order 
to further investigate vdW functionals, we used the so- 
called vdW-optB86b)^i^ which includes the non-local 
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FIG. 3: FeTe energy-volume curves for GGA (open circles) 
and vdW DFT-D2 (full squares) in the PM (upper panel) and 
AFM2 (lower panel) phases. Volume and energies refer to the 
conventional unit cell. Solid lines are the fitted Murnaghan 
curve 



FeSe and FeTe, respectively. 
0.20 
0.15 




FIG. 4: AFMl FeSe (upper panel) and AFM2 FeTe (lower 
panel) c-relaxed total energy GGA (open circles) and DFT- 
D2 (full squares) curves. The solid lines are shown as guides 
for the eyes. 



vdW interaction in the exchange and correlation energy 
functionals. This method was demonstrated to have a 
wide range of applicability and excellent agreement with 
experimental results on different solids in term of lat- 
tice constants, bulk moduli and atomization energies. 
We performed lattice parameters and internal coordi- 
nates optimization varying independently a and c and 
compared the results with DFT-D2 method. The results 
show the complete consistence of the two approaches. 

The use of the corrected c lattice constant has a strong 
effect on band structure and FS near the Fermi energy. In 
order to understand this effect on electronic properties, 
we calculated the band structure and FSs for AFM FeSe 
and FeTe considering both GGA and DFT-D2 relaxed c 
lattice constants. To disentangle the (small) differences 
in the in-plane lattice parameters (see Tab. |l| we calcu- 
lated the equilibrium c constant fixing both a and h ones 
to experimental values. This is a well justified procedure 
to predict interlayer distance widely used in literature in 
the case of layered crystals in which the strongest vdW 
correction comes from the out-of-plane interactionJ^ii^ 

Fig. m show the theoretical results obtained in this 
way for AFM FeSe and FeTe phases which are compatible 
with experiments in a range around 2% and 0.1% for the 



In Fig. [5] and [H we show the band structures calcu- 
lated for the FeSe and FeTe using the above lattice con- 
stants, and in Fig. [7] the relative FSs. Considering the 
GGA c parameter, we observe in FeSe an hole pocket 
at the r point and two very small electron ones along 
the FX direction. In particular, these two pockets are 
related to the presence of a Dirac-like point just below 
the Fermi energy. A nearly Dirac point is also present 
on c*/2 plane, along the ZR line. The electronic states 
changes sensibly using the predicted c lattice constant 
with vdW correction: the hole pocket at the F point is 
now completely filled, while the Dirac point along the FX 
line shifts nearer Ej?, closing all FSs in the FXM plane. 

In FeTe we observe that hole and electron pockets 
along the FX for the GGA c lattice constant transforms 
in electron pocket completely filling the hole one, once 
vdW parameters are considered. 

In conclusion, we studied the effect of the vdW cor- 
rection on the calculation of lattice constants and bulk 
modulus of FeCh superconductors. We showed that the 
vdW correction is fundamental in order to predict lattice 
structure and bulk moduli in agreement with the experi- 
ments, having a large effect on the out-of-plane bonding 
between the Ch-Fe-Ch layers. 
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FIG. 5: FeSe band structure calculated using the GGA 
(dashed lines) and DFT-D2 (solid lines) lattice constants. 




FIG. 6: FeTe band structure calculated using the GGA 
(dashed lines) and DFT-D2 (solid lines) lattice constants. 



sure, strain and surface effects on the structural, elec- 
tronic and magnetic properties of these compounds. 
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FIG. 7: FeSe and FeTe Fermi surfaces (upper and lower panel, 
respectively) calculated using GGA (dashed line) and DFT- 
D2 (solid line) lattice constants. 
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